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The Cygnus 0B2 Association contains some of the most luminous OB stars in our Galaxy and 
the brightest of which are also among the most luminous in X-rays. We obtained a Chandra High 
Energy Transmission Grating Spectrometer (HETGS) observation centered on Cyg 0B2 No. 8a, 
' the most luminous X-ray source in the Association. Although our analysis will focus on the X-ray 

properties of Cyg 0B2 No. 8a, we also present limited analyses of three other OB stars (Cyg 
0B2 Nos. 5, 9, and 12). Applying standard diagnostic techniques as used in previous studies of 
^ . early- type stars (e.g., Waldron & Cassinelli 2001), we find that the X-ray properties of Cyg 0B2 

\^ ' No. 8a are very similar to those of other OB stars that have been observed using high-resolution 

. X-ray spectroscopy. From analyses of the He-like ion fir emission lines (Mg xi, Si xiii, S xv, 

' and Ar xvii), we derive radial distances of the He-like line emission sources and find that the 

. higher energy ions have their lines form closer to the stellar surface than those of lower ion states. 

' Also these /ir-inferred radii are found to be consistent with their corresponding X-ray continuum 

. optical depth unity radii. Both of these findings are in agreement with previous 0-star studies, 

and again suggests that anomalously strong shocks or high temperature zones may be present near 
Q^. the base of the wind. The observed X-ray emission line widths {HWHM ~ 1000 kms""'^) are also 

compatible with the observations of other 0-star supergiants. Since Cyg 0B2 No. 8a is similar 
in spectral type to C Pup (the only 0-star which clearly shows asymmetric X-ray emission line 
profiles with large blue-shifts), we expected to see similar emission line characteristics. Contrary 
53 . to other 0-star results, the emission lines of Cyg 0B2 No. 8a show a large range in line centroid 

shifts (~ -800 to +250 kms^^). However, we argue that most of the largest shifts may be 
unreliable, and the resultant range in shifts is much less than those observed in C, Pup. Although 
' there is one exception, the H-like Mg xii line which shows a blue-shift of -550 kms^^, there are 

problems associated with trying to understand the nature of this isolated large blue-shifted line. 
To address the degree of asymmetry in these line profiles, we present Gaussian best-fit line profile 
model spectra from C, Pup to illustrate the expected asymmetry signature in the residuals. 
Comparisons of the Cyg 0B2 No. 8a best-fit line profile residuals with those of C Pup suggest 
that there are no indications of any statistical significant asymmetries in these line profiles. Both 
the line shift characteristics and lack of line asymmetries are very puzzling results. Given the very 
high mass loss rate of Cyg 0B2 No. 8a (approximately five times larger than previous Chandra 
observed O supergiants), the emission lines from this star should display a significant level of line 
asymmetry and blue-shifts as compared to other OB stars. We also discuss the implications of 
our results in light of the fact that Cyg 0B2 No. 8a is a member of a rather tight stellar cluster, 
and shocks could arise at interfaces with the winds of these other stars. 

Subject headings: X-rays: stars — stars: early-type — stars: individual (Cyg 0B2 No. 8a) — stars: 
winds, outflows — stars: mass-loss — line: profiles — 
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1. Introduction 

The Cygnus 0B2 association (Cyg 0B2) is sig- 
nificant in the history of steUar X-ray astron- 
omy in that Einstein Observatory observations 
of Cyg 0B2 were the first to show tliat early- 
type stars are X-ray sources (Harnden et al. 1979). 
CassineUi & Olson (1979) predicted that OB stars 
should be X-ray sources based on the superion- 
ization stages seen in UV spectra, and suggested 
that the X-rays arose from a coronal region at the 
base of a cool wind. Based on their Cyg 0B2 
X-ray observations, Harnden et al. found that the 
coronal plus cool wind model predicted too low 
a flux of soft X-rays owing to the expected at- 
tenuation by the overlying wind, and they sug- 
gested that the X-rays come from a distribution of 
shocks in the stellar wind. Lucy & White (1980) 
developed the first shock model for these stars, 
and the concept of embedded shocks has remained 
the prevailing picture for the origin of OB stellar 
X-rays. The X-ray observations of Cyg 0B2 us- 
ing Einstein (Harden et al. 1979), ROSAT (Wal- 
dron ct al. 1998), and ASCA (Kitamoto & Mukai 
1996) found that the brightest stars in this asso- 
ciation have significantly larger X-ray luminosities 
in comparison with other OB stars. However, the 
Cyg 0B2 stars have larger bolometric luminosi- 
ties, and the X-ray to bolometric luminosity ratio 
is in line with that of other OB stars (i.e., with an 
observed Lx/Lboi ^ 10"''; Berghofer et al. 1997). 
Among the early type stars in Cyg 0B2, star No. 
8a is particularly interesting because it is one of 
the most X-ray luminous 0-stars (L^ = 1.5 x 10'^^ 
ergs s~^), and because its fast wind has a large 
mass loss rate, {M> 10^''' MQ-yr^^) which is al- 
most as large as that of Wolf-Rayet stars. In ad- 
dition, all of the four brightest stars in Cyg 0B2 
( Nos. 5, 8a, 9, & 12) are also know to be vari- 
able, non-thermal radio sources (Abbott, Bieging, 
& Churchwell 1981; Bieging, Abbott, & Church- 
well 1989). 

Since this association contains some of the most 
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luminous stars in our galaxy (Cyg 0B2 Nos. 9 and 
12), it has been the focus of many massive star 
studies. For example, Knodlseder (2000), using 
the 2MASS near infrared survey data, analyzed 
the extent and population properties of Cyg 0B2 
and proposed that it should be reclassified as a 
"young globular cluster" instead of an OB Asso- 
ciation. As he points out, this association is pro- 
viding opportunities to study the upper end of the 
main sequence of stars in a rather dense popula- 
tion. The subject of massive star clusters has re- 
cently become even more important owing to the 
well received idea that massive stars are formed by 
mergers in tight clusters (Bonnell & Bate 2002; 
Bonneh, Vine, & Bate 2004). The core of Cyg 
0B2 might originally have been such a region. 

The rapidly expanding stellar winds from OB 
stars are often considered one of the most unex- 
pected and important discoveries of NASA's early 
space program (Snow & Morton 1977). Over 
the years we have learned a great deal about the 
wind driving mechanisms and wind properties (see 
Lamers & CassineUi 1999). However, it is the na- 
ture of the processes that lead to shocks and X-ray 
emission that is still unclear and is currently the 
center of attention in OB stellar wind research. 
As shown by recent studies (Waldron & CassineUi 
2001; CassineUi et al. 2001; Miller et al. 2002), 
the resolving power of the Chandra High Energy 
Transmission Grating Spectrometer (HETGS) can 
provide the critical information required to help us 
understand these processes. For the first time, as- 
tronomers are observing spectrally resolved X-ray 
line profiles from OB stars. In particular, tight 
line complexes like the forbidden, intercombina- 
tion, and resonance (fir) lines of He-like ions are 
spectrally resolved by the HETGS. There are rela- 
tively few OB stars that will probably be observed 
with Chandra at this high spectral resolution since 
large exposure times are required, so it is impor- 
tant to glean as much information as possible from 
targets such as the Cyg 0B2 stars. 

Before the launch of Chandra there was only in- 
direct evidence for a shocked-wind origin for OB 
stellar X-rays. These include: a) the lack of suffi- 
cient oxygen K-shell absorption in previous lower 
resolution X-ray spectra (CassineUi et al. 1981; 
CassineUi & Swank 1983; Corcoran et al. 1993); b) 
the consistency of observed O VI UV P-Cygni line 
profiles with an embedded X-ray source (MacFar- 
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lane et al. 1993); and c) theoretical wind instabil- 
ity studies (Lucy & White 1980; Owocki, Castor, 
& Rybicki 1988; Cooper 1994; Feldmeier 1995). 
In our Chandra observations of C Ori (Waldron 
& Cassinelh 2001), C Pup (Cassinelh et al. 2001), 
and S Ori (Miller et al. 2002), we found that far 
more precise information regarding the location of 
the X-ray sources could be determined from anal- 
yses of the He-like ion fir lines (e.g., O vii, Ne ix, 
Mg XI, Si XIII, S xv). In particular, the f/i Hne 
ratio has long been used as an electron density di- 
agnostic for solar-like plasmas. However, for the 
case of OB stars, it is the large mean intensity 
of UV/EUV photons, and not the density of elec- 
trons, that determines the relevant population lev- 
els, as discussed by Blumenthal, Drake, & Tucker 
(1972). The f/i line ratio can be used to derive the 
geometrical dilution factor of the stellar radiation 
field, and this allows us to extract the distance be- 
tween the fir line formation region and the pho- 
tospheric source of UV/EUV radiation. We call 
Rfir the /zr-inferrcd radius of the line formation 
region derived from the He-like ions. 

Our analyses indicate that the derived Rfir 
for the 0-stars correspond reasonably well with 
their respective X-ray continuum optical depth 
unity wind radii as summarized by Waldron & 
Cassinelh (2002). This correspondence in radii 
indicates that the X-rays source regions are dis- 
tributed throughout the wind from just above the 
photosphere out to approximately 10 i?*. How- 
ever, this does not mean that the entire wind is 
emitting X-rays, but rather that there are sources 
of hot X-ray emitting material embedded in much 
cooler gas at all levels of the wind. In fact, the X- 
rays must be arising from only a small fraction of 
the matter in the wind. For example, by assuming 
that the total observed X-ray emission is approxi- 
mately equal to the total intrinsic X-ray emission, 
the resultant emission measures arc found to be 
4 to 5 orders of magnitude smaller than the to- 
tal wind emission measure (Cassinelh et al. 1981; 
Kahn et al. 2001). This small fraction of very hot 
gas is probably in the form of shock fragments and 
filaments in the clumpy and turbulent outflows. 

Although the results seem to support the idea 
that 0-star X-rays arise in shock fragments em- 
bedded in stellar winds, two interesting problems 
have emerged. First, the lines of the highest ion 
stages arise deep in the flow and require shock 



jumps, Av, that arc larger than the local wind 
speed (Waldron & Cassinelh 2001). Several sug- 
gestions have been offered: fast ejecta emerge from 
sub-photosphere regions (Feldmeier, Shlosman, & 
Hamann 2002); in-fall of stalled wind material in 
the form of clumps (Howk et al. 2000); or confine- 
ment of anomalously hot gas in magnetic loops 
(Cassinelh & Swank 1983; Waldron & Cassinelh 
2001; Schulz et al. 2003). Second, the observed 
X-ray line profile shapes do not conform to pre- 
launch expectations. As expected from an expand- 
ing distribution of X-ray sources embedded in a 
stellar wind, the observed X-ray line profiles from 
0-stars with massive winds show a large range in 
broadness with HWHM of ~ 200 to 1500 kms'^ 
However, with the exception of one star, the line 
centroids are not Doppler blue-shifted. Such a 
shifting or skewing of line profiles to shorter wave- 
lengths was expected, as predicted by MacFarlanc 
et al. (1991) as an effect of X-ray absorption by 
the wind. The X-rays from the shocks on the back- 
side (red-shifted) part of the wind should be more 
strongly attenuated than the X-rays from the wind 
material on the front-side moving towards the ob- 
server (blue-shifted). The 04f star ( Pup (the 
only early 0-star spectral type observed so far) 
is currently the only star thus far which clearly 
shows the expected blue-ward shifted skewed X- 
ray line profiles, so it is of special interest to see if 
the 05.5 1(f) star Cyg 0B2 No. 8a shows similar 
line properties. 

The lack of Doppler blue-shifted lines in nearly 
all OB stars is without a doubt the most un- 
expected result obtained by Chandra observa- 
tions. Although the simplest explanation for this 
dilemma is to say that all the X-ray emission 
is located at radii where the wind attenuation 
of the receding X-ray emission is negligible, the 
problem is that these radii are found to be much 
larger (order of magnitude) than the /ir-infcrred 
radii (Waldron & Cassinelh 2001). Even for the 
largest observed Rfir, the wind density is still 
large enough such that most of the red-shifted line 
emission should be attenuated, and blue-shifted, 
asymmetric line profiles should still be observed. 
A possible explanation is that the winds are espe- 
cially porous or 'clumped' which would allow sig- 
nificantly more red-shifted X-rays from the back 
side to propagate through the stellar wind. For 
example, the winds may have small filling fac- 
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tors of absorbing material as would be the case 
for highly clumped winds. The idea of highly 
clumped winds in OB stars is gaining support. 
Waldron & CassincUi (2001) presented the first 
attempt to fit HETGS line profiles with a wind 
distributed X-ray source model and concluded 
that the wind absorption had to be significantly 
reduced either by large reductions in the mass 
loss rate or the X-rays are distributed in a highly 
non-symmetric wind (i.e., clumps). Even for the 
only star showing clear blue-shifts and asymmet- 
ric lines (i.e., ( Pup), Kramer, Cohen & Owocki 
(2003) also found that significant reductions in the 
wind absorption were required to explain these 
profiles. Howk et al. (2000) present arguments 
that the X-rays from OB stars are formed in bow 
shocks around clumps in the winds. The fact that 
shocks lead to highly compressed regions that may 
be Rayleigh- Taylor unstable led Feldmeier, Oski- 
nova, & Hamann (2003) to consider the effects 
of absorption in fragmented winds. They show 
that the fragmented nature of a wind allows the 
X-rays to escape from much deeper in the wind 
than would be the case for smooth, un-clumped 
winds. Interestingly, the idea that clumps could 
be important for the formation and transmission 
of X-rays in OB stellar winds can be traced back 
to Lucy & White (1980) and their first alternative 
model to the coronal idea for X-ray formation. 
The Lucy & White (1980) model had fast clumps 
being driven by radiation forces through the am- 
bient wind and this would produce frontal (or for- 
ward facing) bow shocks. The more recent models 
such as Howk ct al. (2000) have the fast winds 
colliding with slower clumps and these produce 
reverse shocks. Reverse shocks also tend to be 
the dominant source of X-rays in one-dimensional 
hydrodynamic models (e.g., Feldmeier 1995). 

Although wind clumping appears to be a pos- 
sible answer to the blue-shift problem, the clumps 
introduce a dilemma with regards the f/i line 
ratio analysis results. Since a highly clumped 
wind implies that the mass loss rates reported 
for these stars are actually overestimates of their 
true values, then the relationship between the fir- 
inferred radii and X-ray continuum optical depth 
unity radii is no longer correct, and we are then 
faced with a new problem of understanding the 
significance of this observed correlation. In ad- 
dition, whatever level of wind clumping is re- 



quired to explain the X-rays must also be consis- 
tent with the observational constraints established 
at other wavelength bands (e.g., UV, IR, & radio). 
Again, historically, it is interesting to note that the 
main argument against the base coronal model of 
Cassinelli & Olson (1979) and Waldron (1984) was 
that the X-ray absorption was too large to see any 
X-rays arising from regions near the photosphere. 

Our primary target is the brightest X-ray 
source, Cyg 0B2 No. 8a, and will be the focus 
of our discussion, but some information about the 
other OB stars will be presented. In Section 2 we 
discuss the data reduction and present the HETG 
spectra of four Cyg 0B2 stars. Our observational 
analyses for obtaining line emission characteristics 
and the distribution of the X-ray emission are pre- 
sented in Section 3. The conclusions are discussed 
in Section 4. 

2. The Chandra Observations 
of Cyg OB2 

2.1. Data Reduction 

Observations were taken with the Chandra X- 
ray Observatory using the ACIS-S CCD's with 
the HETGS in the optical path. The observation 
of the Cyg 0B2 association [Observation Identifi- 
cation Number (ObsID) 2572] began on July 31, 
2002 at I'' 51™ 04^ UT with an exposure time of 
approximately 65.12 ks. The spacecraft aim point 
was chosen to place the primary target, cluster 
member Cyg 0B2 No. 8a, in the center of the 
focal plane. This gives us a sharp zeroth order 
image and the highest energy resolution spectrum 
possible for this star. During this observation, the 
spacecraft aspect angle was chosen to allow re- 
covery of the dispersed spectra of the three other 
stars, Cyg 0B2 Nos. 5, 9, 12. A "true-color" 
image constructed from the data is displayed in 
Figure 1 which shows the typical rainbow colors 
of the High Energy Grating (HEG) and Medium 
Energy Grating (MEG) dispersed spectra (e.g.. 
Figure 8.1 in the Chandra Proposers Observatory 
Guide, CXC 2002). 

The CIAO tool TGDETECT was used to lo- 
cate X-ray sources in the focal plane image. The 
sources found are displayed in the middle frame 
of Figure 1. A comparison with the Digital Sky 
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survey image of the same region immediately al- 
lows the four brightest sources to be identified 
with Cyg 0B2 Nos. 8a, 5, 9, and 12. Table 1 
lists the adopted stellar parameters for these four 
stars. Although a number of other faint stellar 
sources are detected, this paper concentrates on 
what can be learned from the four HETGS dis- 
persed spectra obtained during this observation, 
primarily Cyg 0B2 No. 8a. Since all four stars 
have large ISM absorption column densities (sec 
Table 1) the extracted spectra reveal essentially 
no information for A > 12 A. The observed count 
rates for Cyg 0B2 No. 8a are 0.102 (MEG±1) and 
0.044 (HEGibl). The only other spectrum strong 
enough to allow detailed analyses of a few strong 
lines is Cyg 0B2 No. 9 which has count rates of 
0.022 (MEGibl) and 0.007 (HEG±1). 

The intrinsic energy resolution of the ACIS-S 
CCDs act as an effective cross dispersion, so in- 
dividual photon events in the dispersed spectra 
can be assigned relatively unambiguously to the 
correct source. When extracting and calibrating 
the dispersed spectra for these four bright point 
sources, the following issues needed to be kept in 
mind: 

• Because these stars are off-axis, the High 

Resolution Mirror Assembly (HRMA) effec- 
tive area is somewhat less and the line- 
response function for Cyg 0B2 No. 5, 9, and 
12 is somewhat broader (e.g.. Figure 8.23 of 
the CPOG, CXC 2002). In particular, note 
the direct image of star Cyg 0B2 No. 5. It is 
so far off the optical axis of the HRMA that 
its zeroth order image is considerably broad- 
ened, greatly reducing the resolution of the 
dispersed spectrum (See pp. 41-50 and 187- 
190 of CXC 2002). 

• As our primary target, the direct image of 
Cyg 0B2 No. 8a is near the center of the 
ACIS-I CCD array. Its whole spectrum for 
the HEG and MEG fell on light-sensitive ar- 
eas of ACIS-S. Because of its location, this is 
also the case for Cyg 0B2 No. 5. However, 
for stars Cyg 0B2 Nos. 9 and 12, their posi- 
tions near the edge of the chip array means 
that part of their MEG -1 and HEG +1 
spectra fell outside the light-sensitive area 
of the chip, resulting in lower effective areas 
at longer wavelengths. 



• When dealing with an observation contain- 
ing multiple X-ray sources, such as this one, 
it is theoretically possible for an individual 
photon event to be at the correct focal plane 
position and energy for more than one source 
(Figure 8.28, CXC 2002). If such a spatial- 
spectral "collision" occurs, the identification 
of a single source for that particular pho- 
ton would become ambiguous. This might 
occur when a MEG spectrum of one source 
overlaps a HEG spectrum of another source 
(see the overlap between the spectra of stars 
Cyg 0B2 Nos. 8a and 5 in our observa- 
tion). To assess how this problem might ef- 
fect this observation, we independently con- 
structed photon event lists for the spectra of 
the four brightest X-ray sources in the field: 
Cyg 0B2 Nos. 5, 8a, 9, and 12. We then 
examined the resulting photon event lists, 
checking for any individual photon events 
which occurred in more than one event list. 
No such photons were found, indicating that 
each of the extracted spectra are basically 
free from contamination by the other sources 
in the field. In each case, the observed sim- 
ilarity of the positive and negative orders 
(when both were available) for each source 
also confirmed that no such contamination 
occurred. 

• The direct images of Cyg 0B2 Nos. 5 and 
12 did not fall on chip S3, the usual pointing 
direction. The default chip set in the CIAO 
tool MKGARF had to be reset to include 
chip 82 to prevent the erroneous calculation 
of zero effective area for the portions of the 
HEG-t-1 and MEG-hl order spectra which 
fell on that chip. 

High-resolution spectra derived from HETG- 
ACIS-S naturally have an extremely low back- 
ground because the intrinsic energy discrimination 
of the ACIS-S effectively acts as a cross dispersion. 
This not only allows photon events at a specific 
position along a spectrum to be sorted to the cor- 
rect order, but also causes most background events 
(which have incorrect energies for a given location) 
to be rejected. 
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2.2. HETG Spectra of the Cyg OB2 Stars 



3. Observational Analysis 



To survey the dispersed spectral data for the 
four bright stars, we co-added the positive and 
negative first order spectra of each of the four 
stars (the 2"'' and 3'''' order spectra were found 
to have too few counts to make any significant 
contributions). The HEG and MEG ±1 count 
spectra for Cyg 0B2 No. 8a are shown in Fig- 
ure 2, and Figure 3 shows the MEG ±1 spectrum 
of Cyg 0B2 No. 9. Unhke other 0-star spec- 
tra, the spectra of these two stars are noticeably 
dominated by only two strong lines, the H-like 
Si XIV line and the He-like Si xiii fir lines. The 
source of this difference from other 0-stars is re- 
lated to the greater high energy emission associ- 
ated with these stars, and the much larger ISM ex- 
tinction. In addition, the Cyg 0B2 No. 8a spectra 
reveals several higher energy ionic line emissions 
(e.g., Ar XVII, Ca xix, Fe xx, Fe xxii, Fe xxiii, 
Fe XXIV, & possibility Fe xxv). Although Fe XX, 
Fe XXII, Fe xxiii, and Fe xxiv were seen in our 
HETGS observation of S Ori (Miller et al. 2002), 
none of these ionic species were seen in our other 
0-star supergiant HETGS observations. 

The other two stars, Cyg 0B2 Nos. 5 & 12, 
have much weaker count spectra. However, we can 
provide a comparison of all four stars by display- 
ing their associated MEG±1 flux spectra in Figure 
4 which provides us with a first order approxima- 
tion of the X-ray line fltixes for comparisons among 
these four stars. For each star, a custom Ancillary 
Response File (ARE) and a Redistribution Ma- 
trix File (RMF) were generated using the standard 
CI AO routines, and the flux spectra were obtained 
by multiplying the count spectra by the bin energy 
and dividing by the associated ARF and exposure 
time. Even though line-shape parameters cannot 
be recovered from many lines in these spectra, the 
fact that individual lines are resolved in the dis- 
persed spectra still opens up avenues of investiga- 
tion which are not possible with CCD-resolution 
spectra such as a non-grating ACIS-S observation. 
However, as discussed in Section 2.1, artificial line 
broadening can occur for off-axis sources due to 
variations in the instrumental line spread function. 
For example, some of the lines in the spectrum of 
Cyg 0B2 No. 5 appear much wider than seen in 
the other stars, but this appearance is most likely 
caused by the large off-axis location of the source. 



3.1. Measuring the X-ray Emission Line 
Properties 

Since the X-ray line emission profiles from OB 
stars are very broad (HWHM ~ 1000 kms^^), 
early-type stars are among the few classes of X-ray 
sources for which line shapes can be resolved using 
the HETGS. This capability allows us to study de- 
tailed line formation processes that we suspect arc 
operating in these stars. Although the currently 
accepted scenario is that a distribution of stellar 
wind shocks are responsible for the observed X-ray 
emission, Chandra observational results are rais- 
ing several key issues concerning the exact nature 
of this distribution, including the possibility that 
not all of the X-ray emission arises from the stellar 
wind. 

A key parameter in determining the structure 
of an X-ray line profile from a wind distribution 
of sources is the stellar wind column density scale 
factor which is proportional to M/{vaoR*) (see 
Waldron et al. 1998). Since the mass loss rate 
of Cyg 0B2 No. 8a is about 5 times greater than 
C Pup (the only 0-star showing well defined blue- 
shifted lines), we find that this column density 
scale factor is ~ 3 times larger than for Q Pup. 
Hence, well pronounced blue-shifts accompanied 
by highly asymmetric line profiles are expected. 

3.1.1. Line Fitting Procedure 

We derive X-ray emission line properties using a 
relatively simple line fitting procedure. We assume 
that all emission lines can be represented by Gaus- 
sian line profiles superimposed on an miderlying 
bremsstrahlung continuum. The line cmissivities 
and rest wavelengths are obtained from the As- 
trophysical Plasma Emission Code and Database 
(APEC and APED; Smith & Brickhouse 2000). 
For each line, the free parameters are the cen- 
troid shift velocity (Vs), the HWHM velocity, 
and the line normalization strength (i.e., the line 
emission measure, EMx)- Since we are only con- 
cerned with the total flux within a line, we find 
that our fits are essentially independent of our 
choice of the continuum temperature. For all fits 
we assume a continuum temperature of 10 MK. In 
this paper, we concentrate only on the H-like and 
He-like lines. For each individual line or line com- 
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plex, wc define a wavelength region large enough 
to cover the emission lines and provide a good rep- 
resentation of the underlying continuum. In many 
cases, the defined wavelength region may include 
other weaker lines, and to ensure we obtain the 
best estimate of the flux in the strongest line and 
continuum, wc include any lines which arc likely 
to contribute to the observed emission. When fit- 
ting a region with multiple lines, we always assume 
that they all have the same Vs and HWHM, and 
only the line and continuum EMx will be differ- 
ent. This approach is clearly justified for the He- 
like fir lines since all three of these lines must be 
formed under the same conditions, and the only 
parameter that could be different is the normaliza- 
tion of each line. Although all H-like lines and the 
He-like i lines are doublets, their separations are 
much less than the resolving capabilities of both 
the MEG and HEG, so the lines are represented 
by only one Gaussian. 

We use the standard statistics to determine 
a goodness of fit. The quoted best-fit values for 
Vs and HWHM represent averages of their re- 
spective 90% confidence region range for each 
parameter, and the quoted errors are associated 
with the difference from the average. These errors 
do not include the associated MEG and HEG in- 
strumental wavelength uncertainties (e.g., see Ta- 
ble 8.1, CXC 2002). The line fiux errors are de- 
termined by the total counting statistics within 
the line. Each model line is folded through the 
HETGS instrumental response functions (ARE 
and RMF) as determined for each star, and the 
best fit values are extracted. To maintain consis- 
tency in fitting the MEG and HEG fines, the 
statistics are based on fitting a bin size of 0.01 
A for both the MEG and HEG which is approxi- 
mately equal to the wavelength resolution limit of 
the HEG instrument (0.012 A; CXC 2002). Be- 
cause the MEG and HEG have different line re- 
sponse functions, we analyzed them separately. It 
should be noted that the fits resulting from the 
two instruments arc consistent with one another, 
especially for the lines with the highest signal-to- 
noise ratio (S/N). 

3.1.2. Determining the Degree of Line Asymme- 
try 

We test the degree of line asymmetry by analyz- 
ing the distribution of the x^ residuals. For exam- 



ple, suppose wc have a triangular shaped emission 
line (e.g., MacFarlane et al. 1991) with the peak of 
the triangle located blue-ward from the rest wave- 
length and the line emission steadily drops from 
this peak value to zero towards the red- ward side 
of the line (there is no emission blue-ward of the 
peak). By attempting to fit this line with a model 
Gaussian line profile, we expect two characteris- 
tic features to emerge: 1) the Gaussian fit will 
underestimate the line ccntroid shift, and; 2) the 
residuals blue-ward of the Gaussian best fit line 
centroid shift velocity will be positive (model un- 
derestimates the data), whereas, red- ward from 
this centroid shift velocity, the residuals will be 
negative (model overestimates the data). Hence, 
a clear residual signature is expected if there is 
asymmetry in an observed emission line. 

To test this hypothesis. Figure 5 shows the 
Gaussian best-fit models for the four strongest H- 
like emission lines of C, Pup which is known to 
have large blue-shifts and asymmetric line profiles 
(Cassinelli et al. 2001). The expected residual sig- 
nature is clearly evident in the Fe XVII and O VIII 
lines, and possibility in the Ne X line. The Mg XII 
line appears to have no line asymmetry based on 
our residual analysis. Hence, we will use these 
C, Pup line fit residual signatures as our benchmark 
for establishing the degree of X-ray emission line 
asymmetries in our Cyg 0B2 No. 8a line profiles. 
Note that this approach cannot be applied to re- 
gions of line overlap such as fir line complexes. 

3.1.3. Best-Fit Line Emission Characteristics 

We apply our line fitting procedure to the four 
H-like and four He- like X-ray emission lines of Cyg 
0B2 No. 8a. The results of our fits are given Ta- 
ble 2 which lists the ion, rest wavelength, observed 
line flux, Vs, HWHM, EMx, the X-ray tempera- 
ture (Ti) range where the line emissivity is within 
75% of its expected maximum value, and the re- 
duced x^ • The EMx are derived as explained 
by Kahn et al. (2001), including ISM absorption 
corrections, and their errors are based on the dif- 
ferences associated with the range in emissivities 
as determined by the range in T^. In Table 2, the 
fluxes for the He-like ion fir lines represents the 
total flux from all three lines, but the listed EMx 
corresponds to the EMx of the r line only. One 
has to be careful in their interpretation of EMx 
derived from the i and / lines since the emissivities 
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used to derive these quantities have not been cor- 
rected for either density and/or UV effects. Hence, 
the r line EMx represents the most reasonable 
estimate since it is essentially unaffected by these 
density and UV effects. The individual fluxes from 
the three fir lines are given in Table 3 (see fir 
discussion in Sec. 3.2). 

The MEG and HEG H-Uke X-ray line fits and 
residuals are shown respectively in Figure 6 and 
Figure 7. Correspondingly, the MEG and HEG 
He-like X-ray line fits are shown respectively in 
Figure 8 and Figure 9. Although the line fits indi- 
cate a tendency for blue-shifted lines (see Vs in Ta- 
ble 2), these shifts are smaller than those of C, Pup, 
which have an average blue-shift of ~ —500 km s~^ 
(Cassinelli et al. 2001). Furthermore, these line 
shifts are significantly smaller than expected as 
discussed in Section 3.1.1, and most of these line 
shifts are similar to other 0-star results which also 
show minimal or no blue-shifts. However, for a 
few lines, large blue-shifts are indicated (see Sec. 
3.1.4). The expected line broadness seen in 0-star 
spectra is also present in the four brightest Cyg 
0B2 OB stars. Our detailed analyses of Cyg 0B2 
No. 8a indicates that all H-like and He-like lines 
have HWHM of ~ 530 to 1100 kms'^ (neglect- 
ing upper limit values), with a suggestion, based 
on the HEG fits, of a decrease at higher Tl. This 
behavior was also noted by Waldron & Cassinelli 
(2002) in their analysis of four 0-stars. However, 
we note that even though a large range of T/, is 
represented in our sample of lines, there is only a 
small relative change in HWHM. 

From the resultant (Table 2) and visual in- 
spection of the line fits, we STiggest that our as- 
sumed input Gaussian line profile model provides 
reasonably good fits to the data. However, before 
we can claim that a Gaussian profile is a realis- 
tic representation of the intrinsic line profiles, we 
must also determine if there is any evidence of 
line asymmetries in the observed line profiles. By 
comparing our best-fit residuals with the C, Pup 
residuals (Fig. 5), we do not see any residual sig- 
natures in the Cyg 0B2 No. 8a lines that would 
suggest the presence of line asymmetry. However, 
there are two possible exceptions, the MEG Si XIV 
and S XVI lines, where both lines show a departure 
from a Gaussian in one wavelength bin (see discus- 
sion in Sec. 3.1.4). Although we stated that this 
residual test is not applicable to the He-like fir 



lines, we point out that there are no obvious resid- 
ual patterns that could be construed as evidence 
of line asymmetries with one possible exception. 
The MEG S xv r line shows evidence of model 
deficient counts in the blue wings of the r line, 
but again, this is not seen in the HEG S xv line. 
Hence, we suggest that since the majority of the 
Cyg 0B2 No. 8a X-ray emission lines do not show 
any evidence of line asymmetries, Gaussian model 
line profiles are consistent with the observed line 
profiles. This is clearly contrary to our expecta- 
tions as discussed in Section 3.1.1, and also sug- 
gests that wind distributed X-ray source models 
(e.g., Owocki & Cohen 2001) will have problems 
in trying to fit these line profiles unless the mass 
loss rate is drastically reduced. In fact, even for 
C Pup a reduction in mass loss rate was also de- 
termined to be necessary in order to explain the 
X-ray line profiles with a distributed X-ray source 
model (Kramer et al. 2003). 

3.1.4- Comments on Peculiar Emission Line 

Characteristics 

With regards to line profile shapes and line 
centroid shifts, our Cyg 0B2 No. 8a results are 
not as straightforward as those obtained from the 
other 0-stars. In our earlier studies, we found that 
the X-ray emission line properties for a given star 
could be categorized as either having asymmetric 
profiles with large blue-shifted lines, or symmet- 
ric profiles with essentially no line shifts. For ex- 
ample, we found that all emission lines in ( Pup 
(Cassinelli et al. 2001) were found to have a large 
blue-shift of ~ —500 kms~^ with many lines dis- 
playing asymmetric line profiles (see discussion in 
Sec. 3.1.2), whereas, for 6 Ori (Miller et al. 2002), 
all lines were found to have symmetric line profiles 
with a small range in shifts (i.e., ±150 kms^^). 
For Cyg 0B2 No. 8a, we find a large range in 
the observed line centroid shifts (Vs in Table 2), 
and no evidence for any line asymmetry (see Sec. 
3.1.3). Although it appears that the Cyg 0B2 No. 
8a X-ray emission line properties do not fit either 
category, in the following discussion we present ar- 
guments which suggest that the line properties of 
Cyg 0B2 No. 8a may not be that unusual, and are 
in fact consistent with what we found for all other 
0-stars, except for C Pup, in having lines that are 
symmetric with minimal blue-shifts. 

From our analysis of the Cyg 0B2 No. 8a X- 



8 



ray emission lines, the range in line centroid shifts 
is -832 to +245 kms~^ with no indication of any 
systematic trends in these shifts (e.g., a temper- 
ature dependence). However, since most of these 
large shifts are associated with the lines from the 
higher energy ions (S XV, S XVI, Ar XVll) which 
have low count rates, implying poorly determined 
shifts, we argue that these large shifts may not be 
real. If we ignore the line shifts from these higher 
ionization state lines, the resultant range in cen- 
troid shifts reduces to -550 to -|-50 kms~^ which is 
still unusual as compared to other 0-star results. 
The cxtremum of this range is due entirely to one 
line, the MEG Mg xii line with a blue-shift of - 
550 kms~^, and without this line, the range in line 
shifts reduces to a value consistent with the major- 
ity of 0-stars. In fact, by averaging the MEG and 
HEG centroid shifts of the strongest lines (Mg xi. 
Si XIII, and Si xiv, neglecting Mg Xll), we find a 
mean shift of only -91 kms^^. 

The real problem is trying to understand the 
large blue-shift of the Mg xii line. In particular, 
a major discrepancy is evident when we compare 
the centroid shift of the Mg xii line with that of 
the He-like Si xill line shift which shows a mini- 
mal blue-shift. The discrepancy in these centroid 
shifts is verified by both the MEG and HEG spec- 
tra. At first, the obvious explanation of this shift 
discrepancy is that the Mg XII line forms farther 
out in the wind than the Si xiii lines, but, based 
on our discussion in Section 3.1.1, the large cen- 
troid shift seen in Mg Xll should also be accompa- 
nied by a highly skewed line profile which is not 
observed (see Sec. 3.1.3). Regardless of this prob- 
lem, the most puzzling aspect is that since both of 
these lines have essentially the same ionic abun- 
dance dependence on temperature, and since both 
of these lines reach their maximum emission at 
~ 11 MK, whatever conditions are responsible for 
forming the Mg XII line must also be forming the 
Si XIII lines which implies that both lines should 
have similar line properties. Although both lines 
have symmetric line profiles, the large difference 
in their centroid shifts is difficult to understand. 
We suggest that some of this discrepancy could 
be related to the low S/N of the Mg xii line [~ 5 
(MEG) and 4 (HEG)] as compared to the S/N of 
the Si xiii lines [~ 9 (MEG) and 5 (HEG)], but, 
at this point, we conclude that the cause of the 
peculiar blue-shift of the Mg Xli line remains un- 



clear. The resolution of this issue will require a 
higher S/N observation. 

In addition, the observed MEG Si xiv and 
S XVI line profiles display possible evidence for a 
narrow highly blue-shifted strong emission com- 
ponent, which would suggest a departure from the 
assumed Gaussian input line profile. Although it 
is intriguing that both MEG lines show this fea- 
ture occurring at the same blue-shifted velocity of 
~ —700 kms~^ (see Fig. 6), which is significantly 
larger that the best-fit Gaussian line model pre- 
dictions for these lines (see Table 2), there are two 
issues that challenge the validity of this feature. 
First, the MEG S xvi fine is very weak. Second, 
this large blue-shifted component is not evident in 
either of the corresponding HEG Si xiv and S XVI 
lines, but this could be related to the lower sensi- 
tivity of the HEG as compared to the MEG. The 
HEG cflFcctivc area in this spectral range is smaller 
than that of the MEG by ~ 44% for Si XIV and ~ 
58% for S XVI. However, we do not wish to ignore 
this feature completely since the shift might be 
real, and could be an indication of an interesting 
new phenomenon occurring in early 0-stars (e.g., 
a high velocity mass ejection from some region of 
the stellar wind or possibility from the stellar sur- 
face). 

3.2. Diagnosing the Location of the X-ray 
Emitting Plasma Using He-Like Ions 

For OB stars the //i ratio provides information 
on the radial distance from the star to the He- 
like ion line formation region (Kahn et al. 2001; 
Waldron & Cassinelli 2001). Even though the He- 
like fir lines often overlap as demonstrated in Fig- 
ure 10, line fitting procedures have been successful 
in extracting the individual fluxes for each of the 
three fir lines. From our studies of these He- like 
line formation regions, we have been finding that 
for most OB stars, the high ion stages are formed 
with Rfir close to the star, which could indicate 
the presence of magnetically confined regions on 
the surface (Waldron & Cassinelli 2001). 

It is important to point out that since the f /i 
ratio is sensitive to the assumed photospheric flux 
model, there is a model dependent uncertainty 
in the derived Rfir (e.g., see discussion by Miller 
et al. 2002). This is particularly relevant for the 
Rfirderived from the Ar xvii, S xv, and Si xiii 
He-like f/i ratios which are sensitive to the un- 
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measurable flux short-ward of 912 A. However, 
since this j ji diagnostic has shown an overall con- 
sistent pattern in the radial locations of the He-like 
X-ray emission lines for several 0-stars, we believe 
that this diagnostic is currently the most viable 
procedure for estimating the radial locations of the 
He-like X-ray emission lines, with the understand- 
ing that the results are sensitive to the assumed 
photospheric flux model. 

Following the procedure outlined by Waldron 
& Cassinelli (2001) we have determined the ra- 
dial positions of the He-like ions of Mg xi, Si xiii, 
S XV, and Ar xvii. The MEG and HEG //t ra- 
tios and derived radial positions, Rjir-, are given 
in Table 3. The //i ratio dependence on radius 
is shown in Figure 11 where the MEG and HEG 
observed ranges in j ji and their corresponding 
radial ranges are indicated by darken curved sec- 
tors. Also shown is our standard X-ray continuum 
optical depth unity plot as a function of radius 
which shows that the Rjir do appear consistent 
with other OB stellar results, i.e., the observed 
He-like X-ray emission is predominantly emerging 
from its associated 'effective' X-ray photosphere 
as first suggested by Waldron & Cassinelli (2001). 
The calculations of these X-ray continuum optical 
depth unity radii are discussed in the Appendix. 
As evident in Figure 11, although the mapping of 
Rfir to the X-ray optical depth unity radii is not 
exactly one-to-one, the overall mapping does ap- 
pear to follow the wavelength dependence of the 
optical depth unity radii reasonably well. One no- 
table exception is that the S xv Rfir is signifi- 
cantly smaller than expected. This could imply 
that the associated high temperatures required to 
produce this emission are only produced deep in 
the wind (< 1.1 stellar radii) and, hence, are ac- 
tually located at a depth where the X-ray opti- 
cal depth is slightly > 2. As seen in other O- 
stars, we also sec that the higher ion stages are 
progressively closer to the stellar surface and, in 
particular, Ar xvii, is essentially on the surface 
(not shown in Figure 11). Figure 11 also shows 
the wavelength dependencies of the X-ray contin- 
uum optical depth unity radii for a mass loss rate 
2 times smaller and 2 times larger than the value 
given in Table 1. The best match between Rjir 
and X-ray optical depth unity radii appears to be 
associated with a mass loss rate that is ^ 2 times 
smaller. In a recent study, Hanson (2003) sug- 



gests that the Cyg 0B2 association may be ^ 35% 
closer than originally though which would imply a 
reduction in the radio determined mass loss rate 
(given in Table 1) by a factor of ~ 2. 

3.3. Temperature Diagnostics of the X-ray 
Line Emitting Plasma 

Of particular interest regarding OB stars is 
the temperature of the X-ray emitting plasma. 
There can be several causes of this X-ray emis- 
sion, ranging from shocks embedded in winds, bow 
shock structures either around clumps in winds 
or around unseen companion stars, or magnetic 
structures near the base of the wind. For example, 
determining the relationship between EM and T 
could provide clues on the X-ray formation pro- 
cess. Here we discuss two line ratio diagnostics 
techniques to derive X-ray emitting plasma tem- 
peratures. 

Miller et al. (2002) were the first to apply the 
H-like to He-like line ratio temperature diagnos- 
tic to an early-type star, the late 0-star 5 Ori. 
In general, one might expect problems using this 
diagnostic for X-rays embedded in a dense wind 
of an OB star since a small change in position 
could produce significantly different wind attenu- 
ations. However, it is becoming increasingly clear 
that each X-ray emission line appears to be arising 
from its associated X-ray continuum optical depth 
unity position. Hence, when taking line ratios, the 
attenuation factors are roughly equal and cancel. 
Normally, the He line used in the H/He line ratio 
is the sum of all fir lines. However, to avoid the 
uncertainties associated with the formation of the 
/ and i lines as discussed earlier, we only use the 
He-like r line in our calculation of the H/He line 
ratio. We have used the APED data to determine 
the dependence of the H/He ratio on tempera- 
ture. Table 4 lists the available H-like to He-like 
line ratios and their derived X-ray temperatures, 
Th/he- Notice that both the MEG and HEG de- 
rived temperatures are in very good agreement. 
Since we expect that the two lines associated with 
each H/He ratio are probably formed in different 
wind locations, our interpretation of this tempera- 
ture diagnostic is that it represents the average X- 
ray temperature between these two locations and, 
in fact, all of these Th/he are found to lie within 
the temperature range specified by their respec- 
tive H and He line Tl values. All of these temper- 
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aturcs along with their associated Rfir support 
earlier claims that, in general, the X-ray tempera- 
tures associated with the majority of the observed 
X-ray emission in OB stellar winds are increasing 
inward toward the stellar surface. 

Another line ratio technique used extensively in 
solar studies is the He-like G-ratio [G = {f + i)/r]. 
Waldron & Cassinelli (2001) were the first to ap- 
ply this technique to an early-type star, the O- 
star C Ori, where they found good agreement with 
Tl- The observed ISM corrected G-ratios for our 
He-like lines are given in Table 3, and the associ- 
ated temperatures, Tfir, are determined by com- 
paring our observed G-ratios with G-ratios cal- 
culated from the APED data. Except for the 
HEG Mg XI G-ratio, there appears to be signifi- 
cant problems with this method as compared with 
Tl and Th/hb- One possible explanation is that 
the r line flux may be reduced by strong res- 
onance line scattering which results in a higher 
G-ratio and a smaller temperature (Porquet & 
Dubau 2000). However, there arc relatively large 
discrepancies in G-ratio values quoted throughout 
the literature. For example, the G-ratio temper- 
ature relation used in the analysis by Waldron & 
Cassinelli (2001) is found to be significantly differ- 
ent than those quoted by Porquet & Dubau (2000). 
As these line ratio discrepancies are resolved, the 
usefulness of this technique for OB stars needs to 
be explored in greater detail. 

3.4. The Emission Measure Distribution 

The range of temperatures in the X-ray form- 
ing regions tells us about the nature of the shocks 
involved. Cohen, Cassinelli, & MacFarlane (1997) 
found that for the near main sequence star r Sco, 
EM{T) oc T~^. Such a power law result is also 
found for bow shock models of clump generated X- 
rays from hydrodynamical theory (Moeckel, Cho 
& CassinelU 2002). In the Moeckel et al. bow 
shock model, it was found that a wide range of 
ionization conditions could be present both be- 
cause there can be very hot matter right at the 
peak of the bow and a whole range of cooler ma- 
terial produced at the oblique shock region of the 
bow. The detection of Ar xvii line emission, for 
example, indicates that there is a hotter source of 
gas in Cyg 0B2 No. 8a than in other 0-stars we 
have studied, such as ( Pup. Moeckel et al. (2002) 
showed that the emission measure versus T distri- 



bution behind a bow shock is a power law where 
EM{T) oc T-4/3. The modeling of r Sco by Howk 
et al. (2000) indicated that the wind is strongly 
influenced by bow shocks. Miller (2002) found a 
rough power law dependence for 6 Ori with a slope 
of —2/3 and it was argued that there could be a 
wind collision with a companion star. However, 
based on the EM values listed in Table 2, there 
is no clear indication of any temperature depen- 
dence and, in fact, one could argue that all line 
EM values are essentially the same. 

3.5. Spectral Fit to the HETG/MEG 

Spectrum of Cyg No. 8a 

In general, fltting the overall HETG spectra of 
OB stars is a diflicult task due to the large range 
in temperatures, the distribution of X-ray sources, 
the widely varying degree of stellar wind absorp- 
tion throughout the wind, the extreme line broad- 
ening, and the quenching of the He-like / line. 
For the case of Cyg 0B2 No. 8a, since there is 
no observed soft X-ray emission due to the large 
degree of ISM extinction, a model with fewer X- 
ray source temperatures can be used. Hence, our 
goal here is to find the simplest model capable of 
reproducing the overall spectral characteristics ob- 
served in the MEG spectrum of Cyg 0B2 No. 8a. 
In addition, this model fit provides us with the 
total observed X-ray flux in the HETGS energy 
band width which can be used for comparisons 
with other observed broad-band X-ray results. We 
use the ROSAT (Waldron et al. 1998) and ASCA 
(Kitamoto & Mukai 1996) derived fitting parame- 
ters as starting points for establishing temperature 
and column density estimates. The X-ray emis- 
sion lines and continuum are calculated using the 
MEKAL emissivity model (Mewe, Gronenschild, 
& van den Oord 1985). In addition, to account 
for the large line broadening and peculiar He-like 
fir line behavior, the emissivity model had to be 
modified by considering: 1) all emission lines are 
assumed to be Gaussian in shape and include a 
pseudo "turbulent velocity" component to mimic 
the line broadening (which is reasonable consid- 
ering that all lines are generally symmetric), and; 
2) artificially large X-ray densities are included to 
suppress the He-like forbidden line emission since 
the current model has no provisions for model- 
ing the effects of an ambient UV radiation field. 
The model fits are then obtained by folding the 
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spectral model through the HETGS instrumental 
response functions (RMF & ARF). The free pa- 
rameters are the temperatures, wind column den- 
sities, and emission measures. The fixed param- 
eters are the ISM column density given in Table 
1, a turbulent velocity of 1150 kms~^, and for 
the low and high temperature components, the 
emissivity densities are respectively 3.2 x 10^^ and 
1.0 X 10^^ cm~^. Note, these densities have no 
physical meaning with regards to the X-ray emit- 
ting plasma since they are only used to simulate 
the correct f /i ratios which are controlled by the 
stellar UV/EUV radiation field. 

We find that a two-temperature, two-wind col- 
umn density model can provide an adequate over- 
all fit to the MEG spectrum. The resultant tem- 
peratures are 3.98 ± 1.67 and 13.88 ± 1.80 MK. 
Their respective wind column densities are (1.22± 
0.60) X 10^2 and (1.41 ± 0.45) x 10^^ cm'^, and 
their respective intrinsic emission measures are 
(3.65 ±2.41) X 10" and (1.06 ±0.24) x lO" cm-^. 
The reduced is 1-29 {x^/DOF = 605/470). 
Comparing these emission measures with the ISM 
corrected line emission measures at similar tem- 
peratures (see Table 2) indicate that these emis- 
sion measures are approximately 17 and 4 times 
larger, and the difference is related to the fact 
that these model derived emission measures repre- 
sent their intrinsic values (i.e., the EMx in Table 
2 have not been corrected for wind absorption). 
These two temperatures and wind column densi- 
ties are consistent with the ROSAT and ASCA 
derived parameters, as well as the observed log 

L^/Lboi of 6.75. Although the ROSAT and 

ASCA fits only required one column density, we 
found that two wind column densities were re- 
quired in order to fit the relative strengths of 
the Mg XII line with respect to the Mg xi and 
Si XIII lines. The majority of the Mg XI emission 
is primarily associated with the 3.98 MK com- 
ponent whereas, the majority of the Mg xii and 
Si XIII line emissions are associated with the 13.88 
MK component. Furthermore, these best fit col- 
umn densities can be used to estimate the loca- 
tion of the two temperature components assuming 
a spherically symmetric wind (sec Waldron et al. 
1998). The resultant radius for the 3.98 MK com- 
ponent is 6.44 ± 2.97 i?*, and for the 13.88 MK 
component, the radius is 4.79 ± 1.38 i?*. Although 
the associated errors in these radial positions in- 



dicate that these locations are essentially indis- 
tinguishable, these model best-fit radial locations 
do suggest that the hotter component is located 
slightly deeper in the wind than the cooler compo- 
nent which is consistent with our fir analysis (see 
Sec. 3.2). In addition, the cool component loca- 
tion is found to fall within the combined MEG and 
HEG Mg XI Rfir range of 2.7 to 6.1 i?*(see Table 
3), and is also consistent with the single compo- 
nent model fit to _R0S'j4Tdata which predicted an 
X-ray location of ~ 4.7 i?, (Waldron et al. 1998). 
The hot component location is found to be larger 
than the combined MEG and HEG Si xiii 
range of 1.7 to 2.4 i?*. The most likely explana- 
tion of this discrepancy is the simplicity of assum- 
ing a two-component fit model. For example, by 
considering a model with a distribution of X-ray 
temperatures around the current hot component 
temperature, significantly more emission over a 
broader spectral energy range would be produced. 
This excess emission would have to be balanced by 
an increase in the absorbing wind column density, 
which in turn would predict a smaller radial lo- 
cation range for this distribution of temperatures. 
In principle this same argument should also be ap- 
plicable to the cooler component. However, since 
most of this soft X-ray emission is masked by the 
large ISM absorption, and this observed soft emis- 
sion has very weak lines with essential no contin- 
uum, it is understandable why a single tempera- 
ture component provides a reasonable fit to the 
softer spectral region. 

The comparison of our model fit with the ob- 
served spectrum is shown in Figure 12 which only 
covers the wavelength region of the strongest lines. 
This shows that the line strengths, broadness, and 
f /i ratios are in very good agreement. The only 
real discrepancy is in the region between 7 and 8 
A where the model has a problem in fitting these 
weaker lines. The Chandra observed log Lx/Lhoi 
is found to be ~ —6.75 which is essentially iden- 
tical to the value determined when last observed 
in 1993 (within - 20 %). Waldron et al. (1998) 
studied the long term X-ray variability of these 
Cyg 0B2 stars and noticed that, of the four, the 
observed X-ray emission from Cyg 0B2 No. 8a 
had remained essentially constant for ~ 15 years. 
We can now suggest that Cyg 0B2 No. 8a X- 
ray emission has now remained constant for the 
last 24 years. In addition, since the Lx/Lboi ra- 
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tio is consistent with the general 0-star behavior, 
we can argue that the X-ray emission processes in 
Cyg 0B2 No. 8a are probably no different than 
those occurring in isolated 0-stars. In particular, 
Chlebowski (1989) found that the Lx/Lboi ratio 
for close binary stars have significantly larger val- 
ues than the observed Lx/Lboi value for Cyg 0B2 
No. 8a. 

3.6. Analysis of the Cyg OB2 No. 9 MEG 

Spectrum 

As summarized by Waldron et al. (1998), Cyg 
0B2 No. 9 is by far one of the most interesting 
variable stellar radio sources where essentially ev- 
ery time a radio observation is obtained it displays 
a different structure. For example, it has displayed 
both thermal and non-thermal characteristics at 
different epochs. However, even when it does ap- 
pear in a thermal state, its radio emission is still 
different, which, assuming frcc-frcc emission, in- 
dicates a highly variable mass loss rate. On the 
other hand, with regards to X-ray emission, Cyg 
0B2 No. 9 is similar to Cyg 0B2 No. 8a in that 
the X-ray emission has remained relatively stable 
for ^ 15 years. In addition, Cyg 0B2 No. 9 is 
also the weakest of the four in terms of X-ray emis- 
sion and the strongest radio emitter, whereas, Cyg 
0B2 No. 8a is the strongest X-ray source and the 
weakest radio source. 

Although Cyg 0B2 No. 9 is the weakest of 
the four main stellar X-ray sources in Cyg 0B2, 
it is the closest of the remaining three to the aim 
point (Cyg 0B2 No. 8a), resulting in the second 
best dispersed spectrum observed. However, the 
HETG spectrum of Cyg 0B2 No. 9 is weak, so de- 
tailed analyses like those for Cyg 0B2 No. 8a are 
not possible. Also, the bigger instrumental line 
spread function makes the emission lines of Cyg 
0B2 No. 9 more difficult to measure. Nevertheless 
by using an RMF-based analysis to study these 
weak emission lines, reasonable reconstructions of 
line shapes can be obtained. The MEG Si xiii and 
Si XIV lines are strong enough to extract informa- 
tion on line profile parameters and line fluxes, and 
reasonable line flux estimates can also be extracted 
for the Mg XI (r line only) and Mg XII lines. The 
following flux values have units of ergs cm~^ s"""^. 
For Si XIV we find a flux = 0.17 ± 0.03 x 10"!^ 
Vs = -100 ± 270 kms-S and a HWHM < 940 
kms~^. For Si xiii we find a total fir flux= 



0.48 ±0.1 X 10-i-\ Vs -700 ±250 kms-^, and 
a HWHM = 620 ± 430 kms"^. The individual 
fir line fluxes in units of 10~^^ are respectively 
0.12 ±0.02, 0.19 ±0.03, and 0.16 ±0.3, which yield 
an //i ratio of 0.60 ± 0.16. Correspondingly, the 
Mg XI r line flux = 0.11 ± 0.02 x 10"^^, and the 
Mg XII flux = 0.08 ± 0.02 X 10-13. 

From this limited amount of information, we 
find three interesting results. First, the H-like to 
He-like line ratios for Mg and Si (see Table 4) pro- 
vide temperatures that are in very good agreement 
with the Th/hc derived for Cyg 0B2 No. 8a. Sec- 
ond, the Si XIII f /i ratio suggests a radial location 
range between 1.4 to 1.9 stellar radii. As discussed 
above, Cyg 0B2 No. 9 had been observed twice 
when it had a thermal radio spectrum. The associ- 
ated mass loss rate estimates from these two times 
are 12.7 x 10"^ Mgyr"! (1983) and 40.0 x 10"^ 
Afgyr-i (1993). Using the 1983 mass loss rate 
value for determining the X-ray continuum opti- 
cal depth unity radius for Si XIII it is found to 
be ~ 1.5 stellar radii. The 1993 value predicts 
an X-ray continuum optical depth unity radii of 
~ 4 stellar radii. Since the 1983 value provides a 
consistent radius with the f/i radius we suggest 
that the 1983 determined mass loss rate is a bet- 
ter estimate of the mass loss rate for Cyg 0B2 
No. 9 at the time of our observations. Third, the 
fit to the Si XIII lines suggests a rather large blue- 
shift of ~ —700 kms~^. If correct the implications 
are very interesting, for we know from the f/i ra- 
tio that the radial location of Si xiii is between 
1.4 to 1.9 stellar radii. These radii correspond 
to a wind velocity range of 915 to 1200 kms~^ 
(using the standard velocity law with a /3 = 0.8; 
Groenewegcn, Lamcrs, & Pauldrach 1989). By as- 
suming a saw-tooth wind shock structure as pro- 
posed by Lucy & White (1980), this wind veloc- 
ity range implies a shock velocity range of 215 to 
510 kms^^ and a corresponding shock tempera- 
ture range of only 0.6 to 3.6 MK, well below the ex- 
pected temperature needed for formation of Si xiii 
(~ 11 MK). We also see a similar situation from 
the Si XIV line and S xvi lines observed in Cyg 
0B2 No. 8a (see Table 2 and discussion in Sec. 
3.1.4). Thus it is difficult to understand these ob- 
servations in the context of this relatively simple 
shock model. 

The line-spread functions of the other two stars 
(Cyg 0B2 Nos. 5 & 12) are too large to be useful 
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for extracting line shape information and fir anal- 
yses which are the focus of this paper. However, 
their MEG flux spectra are included in Figure 4 
which allows one to extract reasonable line flux es- 
timates for the various observed emission lines and 
provides a comparison of the energy flux spectra 
of all four stars. 

4. Conclusions: Are the X- 
ray Properties of Cluster 
Stars Similar to Those of 
Normal OB Stars? 

There are several reasons why the Cyg 0B2 region 
is of special interest. The strongest X-ray source 
in this region, Cyg 0B2 No. 8a, is the second 
example of a very luminous early 0-star and it is 
one of the very few that Chandra will be able to 
observe within a reasonable exposure time. Since 
we are observing a cluster, our one observation 
provided high resolution spectra of three early O- 
stars and one B supergiant. These objects could 
be of special interest because Cyg 0B2 might be a 
young globular cluster. Thus it might be like the 
X-ray emitting regions observed at lower spectral 
resolution in other galaxies. 

The Cyg 0B2 stars are in a relatively tight clus- 
ter that has been called a young globular cluster 
(Knodlseder 2000; Hanson 2003). Hence, it is pos- 
sible, maybe even likely that there is a contribu- 
tion to the X-ray flux from all of the stars from 
wind-wind or wind-star collisions. For example, 
Bonnell & Bate (2002) found that star formation 
simulations indicate that massive stars are gen- 
erally in binary systems which can be relatively 
wide. Thus while it is likely that > 1/2 of all OB 
stars are in binary systems, these could be too 
widely separated for a predominance of X-rays to 
be arising from wind-wind collisions. Now if Cyg 
0B2 No. 8a were to have its X-rays forming from 
such a collision, then we would see evidence for 
this in our wavelength versus radius plots. That 
is, instead of having the radii of line formation lo- 
cated near X-ray optical depth unity, these radii 
would be much larger, comparable to that of the 
separation of early-type stars. A hint of such a 
departure was seen in our study of S Ori (Miller 
et al. 2002), where Si XIII was found to be formed 



well beyond the X-ray optical depth unity radius. 
However, for Cyg 0B2 No. 8a, since the lines 
do seem to be formed near X-ray optical depth 
unity, we do not see a need for postulating any- 
thing other than this star being similar to other 
early-type single stars. 

All previous X-ray observations {Einstein, 
ROSAT, k ASCA) of the brightest Cyg 0B2 As- 
sociation OB stars indicate that these stars show 
clear evidence that the X-ray emission is signifi- 
cantly hotter than most other OB stars. This is 
also evident from our Chandra observations where 
we clearly see Ar xvii and other higher energy 
ions. From this feature along with the fact that 
these stars are very massive, many have suggested 
that these stars are precursors to Wolf-Rayet star 
formation (e.g., Schaller et al. 1992). Unfortu- 
nately, due to the very large ISM extinction for 
the Cyg 0B2 stars, it is very hard to distinguish 
whether the ISM or wind is the dominant ab- 
sorbing component for the soft X-ray emission as 
also noted by Waldron et al. (1998). Based on 
our two component fit to the Cyg 0B2 No. 8a 
MEG spectrum, we find the wind column densi- 
ties required to fit the spectrum need to be larger 
than the ISM column density, suggesting that the 
wind column density is dominant. However, if 
these are precursors to WR stars, then why do 
they have similar standard 0-star behavior (e.g., 
comparable L^/Lboi and consistency with X-ray 
optical depth unity)? Ignace, Oskinova, & Brown 
(2003) point out that X-ray optical depth unity 
radii are expected to be in a range of a 100 to 
1000 stellar radii for WR stars. One main differ- 
ence between an O supergiant and a WR star is 
the stellar radius, and as discussed in Section 3.1, 
the wind column density scale factor is inversely 
proportional to the stellar radius. For example, 
assuming that the mass loss rate and wind speed 
remains the same and the stellar radius of Cyg 
0B2 No. 8a is reduced by a factor of 3 (from 30 
to 10), we find that the location of X-ray optical 
depth unity for 15 A changes from ^ 10 to 50 
stellar radii. Furthermore, we know that WR star 
mass loss rates are typically larger than those of 
the Cyg 0B2 stars which would raise this X-ray 
optical depth unity radius up to several hundred 
radii, consistent with WR results. 

Although these high X-ray temperatures in- 
ferred from our Chandra observation of the Cyg 
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0B2 stars arc unusual for O supcrgiants, they arc 
not overly unusual with regards to all 0-stars. For 
example, the Chandra HETGS observations of the 
main sequence 05 star, 9^ Ori C, shows evidence 
for even higher temperatures and is believed to be 
related to magnetic activity (Schulz et al. 2003). 
In addition, recent XMM-Newton observations of 
the late O supergiant ( Ori (Mewe, private com- 
munication) shows evidence for X-ray emission 
from Ar xvii which was not seen in the HETGS 
observations of this star (Waldron & Cassinelli 
2001). This is interesting since it may be indi- 
cating the presence of a variable high X-ray en- 
ergy component, possibility associated with sur- 
face activity. A ROSAT time series observation 
of C Ori also indicated the possibility that some 
sort of 'flare-like' event was observed (Berghofer 
& Schmitt 1994). 

To summarize with regards to our detailed anal- 
ysis of Cyg 0B2 No. 8a, we do not see evidence 
for the strong wind absorption like that seen in 
WR stars. We also do not see any strong evidence 
supporting wind- wind interactions as evident from 
our fir analysis and the derived Lx/Lhoi- Thus we 
suggest at this point to treat this star as if it is a 
single 0-star X-ray source. Using the single star 
approach we have come to these interesting con- 
clusions. 

a) Two methods were used to derive the radius 
of X-ray formation, the f /i ratio of He- like ions, 
and a fitting of the entire spectrum with a two 
component model with temperatures and column 
densities for each source as adjustable parameters. 
Both methods arc found to predict that the overall 
radial range of the Mg xi and Si xiii line forma- 
tion region is 1.7 to 6.1 stellar radii. 

b) Although this radial range is far below the 
range expected for X-rays from WR stars, this 
does not rule out the possibility that this star is a 
proto-WR star. 

c) X-rays are forming deep in the wind, as was 
the case for ( Pup and C Ori and essentially all 
of the Chandra observed OB stars as shown in 
the temperature versus radius figure of Cassinelli, 
Waldron, & Miller (2003). This demonstrates that 
the simple idea that the X-rays arise from wind 
shocks of about 1/2 the local wind speed is not 
correct. The cause of the hot gas near the stellar 
surface thus remains a puzzle. 



d) Because the Cyg 0B2 No 8a and ( Pup are 
similar early 0-stars with strong winds, we had 
expected that Cyg 0B2 No. 8a would be more 
like C Pup in showing skewed and blue-shifted line 
profiles, but that is not the case. Thus ( Pup still 
appears to be unique among 0-stars. However, 
it also means that there is still a problem in un- 
derstanding of how OB stars can have symmetric 
X-ray line profiles. There are several explanations: 
1) the winds are so porous owing to fragmentation 
(e.g., Feldmeier et al. 2003) that one can see far- 
ther into the wind, even to the far side of the wind 
where red-shifting of line occurs; 2) the X-rays 
are arising from outflowing and in-falling clumps 
(Howk et al. 2000), but their model should only 
be valid for stars with very weak winds that allow 
for stalling and inflow; 3) Ignace & Gayley (2002) 
suggest that resonance line scattering may reduce 
the degree of line asymmetry but their use of a 
Sobolev analysis and its dependence on a smooth 
velocity gradient seems to be a questionable treat- 
ment of an X-ray shock region, and; 4) the X-rays 
are forming in extended magnetic loops, as had 
originally been envisioned by Underbill (1980) and 
we see both hot up flows and down flows in mag- 
netic tubes. 

e) We were able to get X-ray information re- 
garding four 0-stars with one observation that 
had a carefully planned satellite orientation. Al- 
though we were only able to get the most detailed 
information from the primary source, Cyg 0B2 
No. 8a, we demonstrated that we were also able 
to get some useful emission line information from 
one other source, Cyg 0B2 No. 9, and find some 
reasonable line flux estimates for the two weakest 
sources (Cyg 0B2 Nos. 5 & 12). 

Because the OB stars considered here are in 
a tight stellar cluster, we have had to consider 
whether the X-rays were dominated by collisions 
that produce bow shock structures or by processes 
involving enhance magnetic confinement of hot 
plasma. The arguments that we had to consider 
should be of interest to X-ray astronomers study- 
ing galaxies and star burst regions, since they are 
surely looking at clusters of OB stars somewhat 
analogous to Cyg 0B2. In our analysis we have 
developed criteria for deciding on the nature of the 
X-ray emission. In the case of Cyg 0B2 No. 8a, we 
found no evidence that its X-ray emission is dom- 
inated by processes other than those occurring in 



15 



isolated OB stars. 
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A. Calculation of X-ray Continuum Optical Depth Unity 

Analyses of the Ho-likc ion //z lino ratios in OB stars suggest that the derived /zr-inferred radii are 
highly correlated with their respective X-ray continuum optical depth unity radii (i.e., the radial wind 
location where the X-ray continuum optical depth, = 1, for a given wavelength). In general, the stellar 
wind X-ray continuum optical depth for a given wavelength. A, measured outward from radius r is defined 
as 



poo 

(r) = / axir')nH{r')dr' (Al) 

J r 

where ax{r) is the radial and wavelength dependent wind absorption cross section (cm^), and Hh is the 
hydrogen number density (cm^^) defined in terms of the wind mass density as UHir) = p{r)/ jiRmH (where 
vriH is the atomic weight of hydrogen and /i/f is the mean atomic weight). Although we follow the basic 
procedure for determining cr\{r) as discussed by Waldron (1984), several modifications have been incorpo- 
rated. In particular, we have updated the photoionization cross sections (including all K-shell cross sections 
from all ion stages), the coUisional and recombination rates (including dielcctronic recombination), and all 
other relevant atomic data (e.g., elemental abundances) that are available in the most recent version of the 
Raymond & Smith (1979) emissivity code as summarized by Raymond (1988). In addition, the calculation 
of CTA(r) now includes the contributions from 13 elements (H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni) 
and all their stages of ionization. 

To obtain an analytic relationship between radius and ta, two basic assumptions are required. First, it 
is necessary to simplify the integration in eq. (Al). This is accomplished by assuming that (T\{r) can be 
represented as a density weighted average cross section throughout the wind which is only dependent on the 
wavelength and independent of the radial position. This density weighted cross section is defined as 



avoW = TT^. (A2) 



Tx{R 

where tx{R^) and Nw{R^,) are respectively the total X-ray continuum optical depth and wind column density 

(cm~^) as measured from the stellar surface. A representative energy dependent distribution of crnyg[X) for 
the Cyg 0B2 O- stars is shown in Waldron et al. (1998). Hence, tx can be approximated as 

POO 

(A) / nH{r')dr' (A3) 

J r 

where Nw{r) represents the wind column density measure outward from the radial position r. The error 
introduced by using the approximate tx (eq. A3) as opposed to the actual tx (eq. Al) is found to be only a 
few percent. The second assumption requires specification of the wind geometry and wind velocity law. We 
adopt a spherically symmetric wind with a wind speed that is determined by the standard j3 velocity law 
[i.e., v{r) = Woo(l — R*/r)^]. Although this velocity law allows us to obtain a simple analytic expression, it 
produces a singularity in the wind density at the stellar surface. In the actual calculations used to determine 
o'avgW and produce the plots shown in Figure 11, the velocity law is modified by adding a small but finite 
constant velocity term (Vo) to the standard velocity law, and N\y has to be determined numerically. The 
addition of this constant term places an upper limit to the wind density at i?^ , and the value of Vo is typically 
chosen to be equal to approximately one- half of the photospheric sound speed. We realize that for the special 
case of /? = 1, an analytic expression is still obtainable with the addition of Vo, but here we will adopt the 
more appropriate value oi = 0.8 (Groenewegen et al. 1989). With these assumptions, the value of Nw{r) 
is given by (for /3 ^ 1) 

^'^^'■^ " _^°°nff(r-')'^r-' = ^[1 - wiry/^^-'^] (A4) 
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where 



No 



M 



cm 



(A5) 



and w{r) = v{r)/v^. Although our expression for Nyylr) given by cq. (A4) is vahd for cssentiaUy the whole 
wind structure, it will begin to overestimate A''^^ when v{r) becomes comparable to Vq which occurs at 
radii below ~ 1.01-R*. Below this radius limit the overestimate in Nyy is initially a few percent, reaching a 
maximum of ~ 25% at i?*. However, the details in this region are immaterial to our discussion concerning 
the location of X- ray optical depth unity since any depth indicating a location < l.Oli?* can be interpreted 
as occurring essentially on the stellar surface. For additional details, and the case for /? = 1 see Waldron 
et al. (1998). 

Since our goal is to find the relationship between radial position and wavelength for a given value of t\ 
(e.g.jTx), we proceed by first finding an expression for w{r). By combining the results of eqs. (A3) and (A4) 
we have 



w{r) 



I - 



(l-/3)rA(r) 



/3/(l-/3) 



(A6) 



Using the adopted velocity law to obtain the relation between r and w, our desired relationship between 
radius and wavelength for a given t\ is given by 



R* 



(1 - /3)ta 



l/(/3-l)^ 



(A7) 
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Table 1: Adopted Stellar Parameters ^ 




Cyg 0B2 No. 9 


Cyg 0B2 No. 8a 


Cyg 0B2 No. 5 


Cyg 0B2 No. 12 


Spectral Type 


05f 


05.5 


06f+07f 


B8 la 


Tesi (K) 


44700 


38500*^ 


39800 


11200 


log Lfcoi (ergs/s) 


40.19 


39.77'' 


40.01 


39.79 


M (Mq) 


160 


90.5'' 


118 


71 


R{Rq) 


34 


27.9'' 


34 


338 


t)oo(kms~^) 


2200 


2650'' 


2200 


1400 


d (kpc) 


1.82 


1.82 


1.82 


1.82 


log iV/sM(cin-2)) 


22.06<i 


21.92^1 


22.02'' 


22.23'' 


usini(kms~^) 


145 


95'' 


180 


75 


M(M0 yr-i) 


12.7 X 10"*^ 


13.5 X 10"*^ '' 


34.5 X 10-6 " 


38.5 X lO-'^ " 



"Table values from Bieging et al. (1989), unless marked otherwise 
''Table values from Herrero et al. (2002) 

■^Table values from Waldron et al. (1998) - Adopted mass loss rate for Cyg OB2 No. 9 based on our X-ray analysis discussed in 
Section 3.6 

'^Interstellar column density estimates derived from the observed E(B-V) values of Abbott et al. (1984) using the relation 
Nism/E{B -V) = 5.2 X 10^1 cm-2 (ShuU & Van Steenberg 1985) 
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Table 2: X-ray Emission Line Properties for Cyg 0B2 No. 8a 



Ion 


Ao 


line flux/lQ-^^ 
erg cm~^ s~^ 




HWHM 


EM^ 






X'/DOF 




A 


kms~^ 


kms~^ 


10^^ cm 


-3 


(MK) 




Medium Energy Grating 














Ne X 


12.132 


0.27 ±0.04 


0±450 


651 ± 468 


1.36 ±0, 


,50 


6.82 ±2.09 


37.3/46 


Mg xi(/ir) 


9.169 


0.64 ±0.07 


-191 ±380 


998 ± 182 


2.13 ±0, 


.68 


6.71 ± 10.70 


41.0/42 


Mg XII 


8.419 


0.44 ±0.04 


-550 ±314 


1025 ±421 


1.66 ±0, 


.49 


11.45 ±3.51 


16.2/19 


Si xm{fir) 


6.648 


2.20 ±0.16 


-87 ±151 


884 ± 96 


2.90 ±0, 


.88 


11.03 ±3.09 


52.1/40 


Si XIV 


6.180 


1.08 ±0.07 


-36 ±216 


1096 ±316 


2.56 ±0, 


,70 


18.85 ± 6.26 


17.7/14 


S xv{fir) 


5.039 


1.37 ±0.27 


-195 ±640 


< 830 


2.23 ± 1, 


,02 


16.18 ±4.96 


28.3/24 


S XVI 


4.727 


0.36 ±0.08 


245 ± 307 


986 ± 380 


1.42 ±0, 


,55 


29.32 ± 10.49 


10.8/10 


Ar XVll(/zr) 


3.949 


0.30 ±0.10 


-100 ±250 


< 150 


1.56 ± 1, 


.14 


23.74 ± 7.89 


24.6/22 


High Energy Grating 
















Ne X 


12.132 


0.15 ±0.05 


-300 ± 300 


< 1400 


0.76 ± 0, 


,44 


6.82 ±2.09 


22.2/46 


Mg xi(/ir) 


9.169 


0.91 ±0.16 


-146 ±610 


910 ±389 


3.64 ± 1, 


,35 


6.71 ±1.70 


34.5/42 


Mg XII 


8.419 


0.59 ±0.07 


-432 ± 283 


1019 ±320 


2.24 ±0, 


,73 


11.45 ±3.51 


26.4/19 


Si xiii(/zr) 


6.648 


2.45 ±0.28 


-133 ±479 


1017 ±256 


2.57 ±0, 


,89 


11.03 ±3.09 


44.5/40 


Si XIV 


6.180 


1.11 ± 0.11 


50 ±255 


956 ± 285 


2.64 ±0, 


,80 


18.85 ±6.26 


15.9/14 


S xv{fir) 


5.039 


1.86 ±0.41 


-832 ±368 


533 ± 348 


2.15 ± 1, 


,04 


16.18 ±4.96 


21.0/24 


S XVI 


4.727 


0.52 ±0.12 


-35 ± 635 


< 2050 


2.06 ±0, 


,94 


29.32 ± 10.49 


6.2/10 


Ar xvii(/ir) 


3.949 


0.56 ±0.20 


-736 ± 258 


< 650 


3.74 ±2. 


,06 


23.74 ±7.89 


21.1/22 



"T^ range limits corresponds to values in T for which line emission is within 75% of its most probable value 
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Tablo 3: Ho-likn ion fir T^ino Proportins for Cyg OR2 No. 8a, 



Ion 



r flux 



i flux 



/flux 



G-ratio //i ratio 



Tfir 

(MK) 



R* 



Medium Energy Grating 

Mg XI 0.45 ± 0.04 0.15 ±0.02 0.05 ± 0.01 0.48 ±0.07 0.35 ±0.08 > 11 3.13 ±0.39 

Si XIII 1.17 ±0.07 0.46 ±0.04 0.57 ±0.05 0.87 ±0.10 1.20 ±0.17 <5 1.96 ±0.24 

Sxv 0.60 ±0.10 0.41 ±0.09 0.35 ± 0.08 1.26 ±0.44 0.85 ± 0.33 <5 < 1.11 

Arxvii 0.12 ±0.05 0.18 ±0.06 < 0.03 1.57 ±1.04 < 0.16 < 25 ~ 1.00 

High Energy Grating 

Mg XI 0.52 ±0.07 0.22 ±0.05 0.17 ±0.04 0.76 ±0.16 0.80 ± 0.21 5.04 ±2.71 5.13 ±0.93 

Si XIII 1.04 ±0.11 0.61 ±0.08 0.80 ± 0.09 1.34 ±0.23 1.25 ±0.24 < 3 2.06 ± 0.37 

Sxv 0.58 ±0.13 0.67 ±0.14 0.61 ±0.14 2.20 ±0.87 0.92 ±0.35 <4 < 1.18 

Ar XVII 0.28 ± 0.09 0.20 ± 0.07 0.08 ± 0.05 1.00 ± 0.68 0.42 ± 0.37 < 100 ~ 1.00 



"The /, i, r line fluxes are in units 10 erg cm ^ s ^ 
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Tahle 4: H-like to He-like Line Ratio Temperatures for CWe OR2 Nos. 


8a a,nd 9 


Ions 


HETGS/MEG 




HETGS/HEG 






Line Ratio 


TH/He (MK) 


Line Ratio 


TH,He (MK) 


Cyg No. 8a 










Mg xii/ Mg XI (r) 


1.23 ±0.23 


8.46 ±0.52 


1.40 ± 0.36 


8.78 ±0.77 


Si xiv/ Si xiii(r) 


L01±0.13 


12.96 ±0.58 


1.17±0.23 


13.61 ±0.99 


S xvi/ S xv(r) 


0.60 ±0.24 


16.50 ±2.21 


0.91 ±0.42 


18.90 ±3.32 


Cyg No. 9 










Mg xii/ Mg xi(r) 


0.94 ±0.39 


7.70 ± 1.02 






Si xiv/ Si xiii(r) 


1.18 ±0.42 


13.60 ± 1.74 
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Fig. 1. — This shows a "true-color" X-ray image showing the Chandra spectra of the central region of Cyg 
0B2. For all images, North is up. Photons are color-coded according to the following scheme: Red < 1.5 keV, 
Green = 1.5 — 2.5 keV, and Blue = 2.5 — 8 keV. The middle panel shows the ACIS-I focal plane image with 
the sources (shown by green circles) found using the CIAO tool TGDETECT. The bottom panel shows the 
same regions on a Digital Sky Survey of the same region. Although some of these detected sources correspond 
to stars in the Cyg 0B2 Association, several circled regions shown in the middle and bottom panels are not 
actual X-ray sources. For example, the straight rows of circled regions along the MEG spectrum on either 
side of Cyg 0B2 No. 8a and a few other circled regions (as determined by the automatic source detection 
routine) are simply marking the locations of bright emission lines in the dispersed X-ray spectra of the four 
brightest stars. 
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Fig. 2.— Observed HEG and MEG il'** order count spectra for Cyg 0B2 No. 8a. The most likely line 
identifications are indicated. The bin size is O.Of A. 
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Fig. 3.— Observed MEG ±1"* order count spectrum for Cyg 0B2 No. 9. The most likely line identifications 
are indicated. The bin size is 0.01 A. 
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3 4 5 6 7 8 9 10 11 12 13 
wavelength A 

Fig. 4. - Comparison of the combined MEG ±1"* order flux spectra for all four stars (Cyg 0B2 Nos. 8a, 5, 
9, & 12). The rest wavelengths of the H-like and He- like lines are labeled. The off-axis degradation in the 
flux spectra of Cyg 0B2 Nos. 5, 9, & 12 is clearly visible below 6 A. The bin size is 0.02 A. 
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Fig. 5. — Shows comparisons of the four strongest Zeta Pup H-like best-fit line model spectra (solid- line) 
with their corresponding observed MEG ± 1 count spectra and their associated residuals. The horizontal 
axis is expressed as the velocity shift relative to the rest wavelength of the line. The vertical dashed-lines 
in the residual plots represent the best-fit line centroid shift velocities {Vs). In order to obtain reliable flux 
measurements of the line and continuum emissions for Ne X, the input model spectrum included the weaker 
Fe xvii(~ 12.26 A) located at ~ 3200 kms'^ The bin size is 0.01 A. 
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Fig. 6. — Shows comparisons of the Cyg 0B2 No. 8a H-like best-fit line model spectra (solid-line) with 
their corresponding observed MEG ± 1 count spectra and their associated residuals. The horizontal axis 
is expressed as the velocity shift relative to the rest wavcilciiigth of the line. The horizontal dashed-line 
in the best-fit plot represents the predicted model counts for the continuum. The vertical dashed-lines in 
the residual plots represent the best-fit line centroid shift velocities (Vg). In order to obtain reliable flux 
measurements of the line and continuum emissions for Ne X, the input model spectrum included the weaker 
Fe xvii(~ 12.26 A) located at ~ 3200 kms'^ The bin size is 0.01 A. 
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7. — Same as Figure 6 except for comparisons with the HEG ± 1 count spectra. 
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Fig. 8. — Shows comparisons of the Cyg 0B2 No. 8a He-hke best-fit line model spectra (solid- line) with 
their corresponding observed MEG ± 1 count spectra and their associated residuals. The vertical dashed- 
lines indicated the rest wavelengths of the fir lines. Note, although the i line is actually a doublet, each 
input model spectrum assumes only one line. The bin size is 0.01 A. 



32 




9.1 9.2 9.3 
wavelength A 



6.6 6.7 
wavelength A 



15 



10 



5 - 




2 

OH 



I ' ' ' ' I 
sxv 



-H-4- 



I .... I 



I I I I I I I 



-H 

o 
w 

X 



4.95 5.00 5.05 5.10 5.15 
wavelength A 




5.20 3.85 



3.90 3.95 4.00 4.05 
wavelength A 



4.10 



9. — Same as Figure 8 except for comparisons with the HEG ± 1 count spectra. 
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Fig. 10. — Shows the best-fit input model spectra for two He-likc fir line triads (Mg xi& Si xiii) prior 
to folding through the MEG ± 1 instrumental response functions (ARF & RMF). The input spectra are 
normalized by their respective continuum. The contributions of each line to the total (black line) line emission 
is indicated: r-line (blue) , i-line (green) , and f-line (red) . This demonstrates that the high energy resolution 
capabilities of the HETGS provides a clear distinction of each line's contribution to the overall emission, and 
allows us to extract individual fir line characteristics. 
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Fig. 11 — Top - The Cyg 0B2 No. 8a HEG and MEG He- like f/i ratio dependence on radius. The observed 
range of the f/i ratios and their associated radial ranges are indicated by the broad darkened line sectors. 
Bottom - The corresponding X-ray continuum optical depth unity (solid line) wavelength dependence on 
radius (assuming the mass loss rate given in Table 2) and the observed range in radii associated with each 
He-like f/i ratio. In all plots the Ar xvii radial range is not shown since it is below 1.001. The dashed lines 
correspond to a factor of 2 increase (larger radii) and a factor of 2 decrease (smaller radii) in mass loss rate. 
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Fig. 12. — Comparison of the best-fit model spectrum (red line) to the MEG if* order spectrum for Cyg 
0B2 No. 8a. The spectrum only covers the wavelength region of the strongest lines. The corresponding 
residuals of the fit are shown in the bottom panel. The inset highlights the wavelength region of the two 
strongest lines, Si xiv and Si xiii. The bin size is 0.02 A. 
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